Pressure induced FFLO instability in multi-band superconductors 
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Multi-band systems as intermetallic and heavy fermion compounds have quasi-particles arising 
from different orbitals at their Fermi surface. Since these quasi-particles have different masses or 
densities, there is a natural mismatch of the Fermi wave- vectors associated with different orbitals. 
This makes these materials potential candidates to observe exotic superconducting phases as Sarma 
or FFLO phases, even in the absence of an external magnetic field. The distinct orbitals coexisting 
at the Fermi surface are generally hybridized and their degree of mixing can be controlled by 
external pressure. In this Communication we investigate the existence of an FFLO phase in a 
two-band BCS superconductor controlled by hybridization. At zero temperature, as hybridization 
(pressure) increases we find that the BCS state becomes unstable with respect to an inhomogeneous 
superconducting state characterized by a single wave- vector q. 

PACS numbers: 74.20.-z, 74.25. Dw, 74.62. Fj, 74.81.-g 



Asymmetric superfluidity refers to Cooper pairing in 
systems with mismatched Fermi surfaces. This phe- 
nomenon comprises the FFLO^ type of ground states 
where the mismatch between bands with different spin 
orientations is produced by an external magnetic field. 
It also occurs in cold atom systems where the mismatch 
is due to a different numbers of interacting fermionsSii^. 
Besides, it may appear in the interior of neutron stars 
where the pairing of up and down quarks in different 
numbers can give rise to color superconductivity (see for 
example, Refsi^i^). 

In multi-band metallic systems as inter-metallic com- 
pounds and heavy fermions, electrons arising from 
distinct atomic orbitals coexist at a common Fermi 
surface^'^. Since these electrons have different effective 
masses or occur in different numbers per atom, there is 
a natural mismatch of the Fermi wave-vectors of these 
quasi-particles. As a consequence, we may expect to find 
the physics associated with asymmetric superconductiv- 
ity in these systems even in the absence of an external 
magnetic field. In these materials, the wave functions 
of electrons in different orbitals hybridize due to their 
overlap. In particular, the mismatch of the Fermi wave- 
vectors is affected by hybridization. Since pressure con- 
trols hybridizatior^, we show that in multi-band super- 
conductors it plays a role similar to that of an external 
magnetic field in the study of FFLO phases. It has the 
advantage that this pressure induced FFLO phase does 
not compete with orbital effects which arise when apply- 
ing an external magnetic field to a superconductor. 

The problem of superconductivity in systems with 
overlapping bands was treated originally by Suhl, 
Matthias and Walker—. However, these authors did not 
consider inter-band pairing as this is negligible in the case 
where the critical temperature is much smaller than the 
effective inter-band splitting. 

Recently, we have investigated asymmetric supercon- 
ductivity in multi-band metallic systems in the presence 
of intra and inter-band interactionsiS. We have studied 
the different types of homogeneous ground states which 



appear as hybridization is changed. In the inter-band 
case, as hybridization increases there is a first order tran- 
sition from the BCS state^^ to the normal state. In be- 
tween these states there is a gapless metastable phase 
with similarities to the Sarma phase^^ which has had re- 
newed interest in recent yearS'^*i^. The instability of the 
BCS state is related to the appearance of a soft mode 
at a characteristic wave-vecto r^°'^'^ . This suggests that 
an alternative ground state as hybridization increases is 
an inhomogeneous superconductor of the FFLO type. In 
this paper we investigate the existence of such state. Dif- 
ferently from Ref;^, we consider the situation where the 
dispersion relations of the bands overlap at the Fermi 
surface, such that, their Fermi wave-vectors are equal. 
In this case inter-band interactions must be taken into 
account. 

The effective Hamiltonian describing the two-bands 
metallic system, hybridization and pairing of quasi- 
particles with a net momentum q is given by, 
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where the inhomogeneous superconducting order param- 
eter is. 



(2) 



with g the strength of the attractive interaction. The 
dispersion of the quasi-particles is given by, 
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where, 



(fc) = a^k'^, 
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and a < 1 is the ratio of the effective masses. 

The Green's function method is used to obtain the 
BCS-hke order parameter 



(5) 



where /(w) is the Fermi function. 

In order to calculate the relevant Greens functions we 
obtain their equations of motion. In particular for the 

anomalous Greens function ^^afc_|_2 ; this is 

given by, 
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After some long calculations we obtain for the anoma- 
lous Green's function, 
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and 
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The poles of the Green's function, D{uo) — 0, in Eq. [9] 
yield the excitations of the system. Substituting the dis- 
persion relation of the bands. 
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ak^ + a— ± a fc .~q 
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in Eq[9l we obtain a complete fourth degree equation for 
the energy uj of the excitations, 

(10) 
(11) 



where 
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with vp the Fermi velocity and we have neglected terms 
of 0[q^) as usual. 

In order to solve this equation we introduce the change 
of variable, 



UJ —> u — -^ = u + vpq- 2 



cos 9 
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which yields a depressed equation of the fourth degree 
-I- /3u^ + 7u + A = (13) 



where 
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(7) up to linear terms in q. In the case V — 0, a = 1, 
e% = e\: the fourth order equation reduces to a product 
of two identical second order equations. The roots of this 
second order equation yield the excitations found in the 
usual FFLO problem. 

The problem above is still quite intractable. This is 
due to the different masses (a ^ 1) of the quasi-particles 
that in combination with mixing has a very strong desta- 
bilizing effect on the FFLO state. The effects of hy- 
bridization are stronger at the points in fc-space where 
the bands cross, i.e., for ej! = ■ Analytical progress 
can be done if we assume the case of homotectic bands, 
i.e., we take e\ = ae^, and = €k- The crossing of 
the bands takes place exactly at the Fermi surface, at 
= 0. Furthermore, to make analytical progress we 
consider that the ratio between the masses of the quasi- 
particles a is very close to unity, i.e., we write a = 1 — e, 
and neglect terms of order e^. In this case we can find a 
solution for the depressed fourth order equation given by 
Eq. M 

The energies of the excitations in this case are given 
by a; = cijj^(A;), where, 



(14) 



with, 



Wi2(fc) = V^fei V^fe. (15) 

The quantity J/i = -fe/4 = w_Fg[(l + a) /2] COS0. Also, 

A, ^{1- e)el + V' + Al + Oief (16) 

and 

B, ^W'[{l^e)el + Al] + 0[er. (17) 
These equations yield 

uJi2{k)=^k±V (18) 
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where 



A2. 



When calculating the gap function Ag we find, after a 
change of variables, the following integral, 



1 /■+°° 

Gk{S^i) = — du)Dx{uJ + Sfi) Im 
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where /(w) is the Fermi function, Dx{u>) is given by Eq. 
[8] above and the denominator of the anomalous Greens 
function is given by: 



Using that. 
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UJ — UJ2 LO + UJ2 
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Recalling that in the equation above, lu ~* uj + ie, and 
taking the imaginary part, we obtain that Gk{SfJ.) is a 
sum of three terms, Gk{Sfj.) = G\{5 ^i) + Gl{5 ^jl) + G\{5 ^) 
with, 



Gl{5^i) = ^ 



^|2-X:[/(i?L) + /(£^L)]|, 



Gl{5^^) = 



Y.{-iy-\f{Ei^) + .f{Ei_)] 



and 



j=l,2 

where = Cfe + c^{V + (5^^) and i;^^ = Cfc + (^(V^ - (^A*) 
with a — ±. We have omitted terms of 0{q)^ and 0(e)^- 
When calculating the gap equation A, = Gk{5fi) at 
zero temperature, the Fermi functions are expressed in 
terms of 9 functions and this imposes severe restrictions 
on the sums over k . When these sums are performed 
and angular integrations are carried out, the only contri- 
bution which remains is that arising from G'J,(5/z). The 
gap equation can finally be written as, 
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Subtracting the T = gap equation for a BCS supercon- 
ductor, 
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with awl, from the left hand side of Eq. [191 we obtain 
in the weak coupling approximation. 
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where p is the density of states at the Fermi level. The 
integrals over k (J dk) on the right hand side are per- 
formed taking into account the constraints imposed by 
the 9 functions. They yield. 
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This equation has real solutions only if -I- Vpq > Aq 
where Vp = vp {I + a) /2. Let us consider the case a = 
— 1, Recalling that 5fi = Vpq cos 9, the integral above can 
be rewritten as, 



9P 



Ay/a 2v*pq 



dx sinh 



where we used the change of variables, x — —v*pq cos 9. 
In fact the integrals are independent of a and the result is 
simply twice that for a given sign. Respecting the limits 
of integration in different cases to obtain a real result, 
the final gap equation is given by, 
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where G{x) is the function - , 



(21) 



G{x) = X cosh X 
= 0,\x\ < 1 
= -G{-x),x < 



1, Ixl > 1 



Notice that the mass ratio a cancels out explicitly in the 
gap equation, Eq. [211 It's role at least for a w 1 is just 
to renormalize the Fermi velocity. From this equation 
we find that for the FFLO state to be a solution it is 
necessary that q — q/{V/vp) > 1. Also, since G{\x\ < 
1) = 0, the solution for V < Vf{q) = Ao/{l + q) is always 
Aq — Aq, i.e., the BCS state. Thus a necessary condition 
for the FFLO state is y > Vf{q). The upper critical 
value of the hybridization V2{q) below which the FFLO 
state can be a solution of the gap equation is obtained 
taking the limit of Eq[2Tlfor Aq ^ 0. The results can be 
expressed 
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FIG. 1: Vi (dashed) and Vi as a function of the reduced 
wave- vector. 

In Figure [1] we plot Vf{q) and Vfiq) as a function of 
the reduced wave-vector and it is clear that there is a 
range of values for the hybridization < V < for 
which a FFLO phase may exist. The maximum value of 
V2 occurs for q — q^ « 1.2, which when substituted in 
the equation above yields Vc — V2^(9c) ~ 0.75Ao. This 
value of q is that which minimizes the free energy in 
the range of stability of the FFLO phaseiiii. The value 
Viijl) above marks the limit of stability of the FFLO 
phase. The actual value of the hybridization for which 
the first order phase transition occurs is obtained consid- 
ering the energies of these states. The argument is similar 
to that of Chandrasekhar and Clogston^^ to obtain the 
critical field in BCS superconductors. Here we have to 
consider the hybrid bands. In the limit of very small 
mass differences their dispersion relations can be easily 
obtained and are given by, loi^2 — [(1 + a.)/2]ek ± V. 
On the other hand the condensation energy for a sys- 
tem of unequal masses was obtained in Ref.~. This is 
similar to that of a system of identical particles with 
the mass m replaced by 2mr, where the reduced mass, 
rrir — mafrih/ {nia + ttl},) — ma/{l + a) in our nota- 
tion. The chemical potential is also modified and given 
by, = [iJLa + /i6)/2 = [(1 + a)/2\^ia- Then the effec- 
tive particles have dispersion, — [(1 -I- a)/2]efc. Com- 



paring the condensation energy of these quasi-particles, 
Ec = (1/2)p*Aq with the energy associated with hy- 
bridization, Ev — p*V^, one obtains a critical hybridiza- 
tion, Vc = A0/V2 « O.7IA0, above which BCS super- 
conductivity becomes unstable. In these expressions, 
p* is the density of states at the Fermi level of parti- 
cles with dispersion relation = [(1 -I- a)/2]ek- Conse- 
quently there is a window of values for the hybridization 
(O.7IA0 < V < O.75A0) where we can expect a FFLO 
phase to occur. The transition at V2 is a continuous sec- 
ond order transition from the FFLO to the normal state. 

Our results have a close similarity to the usual FFLO 
approach for a superconductor in an external magnetic 
field. This was anticipated from the form of the dis- 
persion relations, Eqs. 1181 where V enters formally as 
an external magnetic field. However, the analogy with 
the usual FFLO stops there. The Greens functions in 
the present case have four poles, instead of two and the 
numerator of the anomalous Greens function (Eq. [S]) is 
much more complex and includes an angular dependence. 
At the level of the Hamiltonian, Eq. [Tl V mixes different 
states and from this point of view it acts like a transverse 
field and not as a polarizing longitudinal field. The latter 
only repopulates the states while the former changes the 
nature of the quantum states. 

The FFLO phase in condensed matter systems has long 
been sought. Here we point out the possibility of attain- 
ing an inhomogeneous superconducting state by applying 
pressure in a multi-band superconductor. The existence 
of quasi-particles belonging to different orbitals in a com- 
mon Fermi surface provides a natural mismatch. It can 
be controlled by pressure and this, as we have shown, 
offers the possibility of finding new inhomogeneous su- 
perconducting states tuning this external parameter. 
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